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Abstract: 

We analyzed two form factors T\(q 2 ) and T2{q 2 ) which enter the B — » K*£ + £~ decay rate 
and are accessible in the lattice calculations for B — > K*^. The extrapolation to q 2 = 0, 
necessary for B — > K*7 is examined. The flat behavior of T 2 (q 2 ) is shown to be disfavored. 
From the double constrained unitarity bounds by the lattice result at small recoil and the 
light cone sum rule one at large recoil, the pole/double pole fit is preferred with the 'pole' 
masses as free parameters. We also show a consistency of two approaches in the heavy 
quark limit. Besides several byproducts of this analysis, our final (conservative) result for 
the form factor that determines B — > K*^ rate is T(0) = 0.15 ± 0.04. 
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1 Introduction 



In recent years there has been an enormous effort put into studying the flavor changing 
neutral current (FCNC) decays. Unlike the other rare decays, B — > K*^y is more attractive 
because its branching ratio is proportional to G 2 F a and not to G 2 F a 2 . Indeed, experi- 
mentally it was detected and measured. Today, the increasing statistics allows better and 
better accuracy so that the 10% of experimental error seems to be achievable quite soon. 
Last summer, CLEO [Q reported: 

B(B -> K*j) = (4.2 ± 0.8 ± 0.6) x 10~ 5 (1.1) 

which makes around 20% of the inclusive branching ratio. Based on b — > 57, this decay is 
remarkable from several reasons: it allows the access to the V ts Vtb combination of Cabibbo- 
Kobayashi-Maskawa matrix elements; at lowest order in the Standard model (SM) it occurs 
through loops and since the contribution of W boson is small, one expects signs of the non- 
SM physics. That is the main reason why one tries to calculate this decay as precisely as 
possible in the SM. Once the decay rate is fully determined in SM, a discrepancy in its 
comparison with experiment could give the strength of the contribution that comes from 
the physics beyond. The inclusive rate is much more complicated and the results suffer 
from large theoretical uncertainties. One relies then on exclusive process. Theoretically, a 
series of quark models 0, ^ [4[] and phenomenological approaches [|5] were employed. At 
the level of guidance and the phenomenological understanding they are very useful. But, 
the correct estimate of the decay rate should be based on QCD. 
This decay is described by the effective hamiltonian: 

AC 8 

H eff (b -> 57) = -^-VZVt^QiriO, (1.2) 

which is obtained by systematically integrating out heavy degrees of freedom from the 
original theory while descending on the scales of interest (see for example H 0). Oi(n) 
represent the complete set of operators and Cj(/i) corresponding Wilson coefficients at the 
low hadronic scale \i. The dominant contribution to B — > K*^ decay comes from the O-j 
operator (electromagnetic penguin) which is identified as a short distance contribution. The 
other contributions (long distance) were first questioned in || and come from the vector 
dominance possibility B — > K*J/ifj(ip'...) —>■ K*j. In fact, next to leading calculations 
|| show that the contributions of other operators are bigger than it was anticipated at 
leading order. The significant long distance part in the case under study comes from the 
operator O2 ||10|| , when a photon is emitted from the charm lineQ. This part was recently 



calculated by the classical QCD sum rules [11], and estimated to be a small effect for 



the branching ratioQ. The step which is absolutely necessary is to precisely evaluate the 
dominant contribution, and we therefore turn our attention to Oj. It is convenient to 



1 For B — ► p7, long distance effects are far more important. See 
2 The mixing with O2 operator was also studied for the inclusive case 
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parametrize the hadronic matrix element for B — > K*^y as: 
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where q = p — k and Ti^iq 2 ) are the form factors. T 3 (g 2 ) does not contribute B — > K*j 
decay rate and will not be considered in what follows. From 07^75 = — fe^ a/ 3<T a/3 one 
obtains Ti(0) = T2(0). This condition is quite important for lattice analyses since it 
represents the important constraint for the extraction of the form factor at the physical 
point for the on-shell photon. The decay rate is then given by: 
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(1.3) 
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where we wrote Ti^(0) = T(0). Now, it is obvious how important it is to have the exact 
value of the form factor. To implement the QCD in the calculation of T(0), we must have 
the knowledge of the nonperturbative part of the hadronic matrix element. As it is well 
known, only the lattice QCD and the QCD sum rules can give some answer on that part. 



2 Penguins on the lattice 

A huge effort to calculate this hadronic matrix element was done by lattice groups. To 
date, there are four lattice groups [|13|, [14], [15|, [16|] that reported results on the T(0). 
Recall that lattice calculations involve numerical simulations: using the euclidean space 
functional integral formulation of the theory, corresponding correlators are calculated on 
finite lattice (L, a) so that the functional integral becomes multidimensional and can 
be handled by Monte-Carlo methods. This in principle allows us to calculate every- 
thing from the first principles of QCD0. But the problems are multifold: e.g. most 
of the calculations are performed in the 'quenched' approximation in which the inter- 
nal quark loops are left out; the actual lattices cannot accommodate large quark masses 
(a- 1 ~ 2 GeV(/3 = 6.0), 2.8 GeV(/3 = 6.2), 3.6 GeV((3 = 6.4)). We also have to inject 
momenta to B and K* to cover the kinematically accessible q 2 , which requires ever larger 
lattices. Problems related to the treatment of fermions on the lattice are well known. The 
explicitly introduced (Wilson) term to the naively discretized fermionic part of the QCD 
action, introduces discretization errors of order 0(a). Most of the today's simulations 
are performed with the Wilson action and with its improved version i.e. Sheikholeslami- 
Wohlert (SW) action, where the discretization errors are reduced to O(gQa). However, 
these errors can be quite disturbing in the studies of the heavy flavor transitions where 
the corresponding matrix elements of bilinear operators pick up the lattice artifacts of or- 
der 0(ma). Very encouraging is the recent proposal of the ALPHA collaboration for the 



'For a good review, see 17 
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improvement of the SW action (see for example [T^fl), by which the discretization errors 
are reduced to only 0(a 2 ). At the same time, this proposal also allows us to reduce the 
errors for bilinear operators. Originally, the improvement was carried out in the chiral 
limit, which was very recently extended to the case of non-zero quark masses in [Tj|. Con- 
sequently the errors on form factors values out of zero recoil region are expected to be 
considerably reduced. This program should be used in future simulations^. 

Still the problem of the small heavy mass on the lattice is ubiquitous and the results 
obtained in such a situation must be extrapolated to the physical m#. In our case, this 
can be done with the help of heavy quark symmetry which offers important scaling laws: 

HG&J ~ M 1 ' 2 , T 2 (q 2 max ) ~ ^ (2-1) 

The pole dominance prescription then, for the q 2 dependence, gives: 

T ^ = tt 2 ^. < 2 - 2 > 

P 

which is compatible with the above scaling laws and Ti(0) = T 2 (0), if T 2 (q 2 ) is flat (n = 0) 
and Ti(q 2 ) being dominated by a pole (n — 1), or if T 2 (q 2 ) is pole like and Ti(q 2 ) double 
pole (n = 2), and so on. The actual precision of the lattice data does not allow for a 
definite conclusion on this matter. The reason is clear: Directly 'measured' values of form 
factors are obtained for smaller masses and limited range of transfer q = p— k. Even with 
the bigger statistics, a firm conclusion on the behavior (|2.2|) was not claimed. To some 
extent, the shapes ( |2.2|) can be tested to 'interpolate' to q 2 = since the direct results are 
obtained around this point. Again, to make these tests more reliable, lattice data should 
be more precise for several q ^ and it is particularly at this point that the new program 
for the improvement is expected to play an important role. The fitted values of the form 
factors at q 2 = should be then extrapolated from lattice accessible heavy mass to the 
physical tub- The problem was the lack of such a scaling law. In the case of B — > p{js}lv 
this scaling law can be obtained in the heavy quark limit from QCD [^, ^l|] which turns 
out to be ~ M -3 / 2 for all form factors. The situation for B — > K*j is (apart from small 
SU(3) breaking, unimportant for these scaling laws) the same [^| and for future analyses, 
these scaling laws must be used. As a cross check of the analysis, the form factors are 
first extrapolated to and then to q 2 = by means of ( |2.2| ) and with the help of the 
constraint Ti(0) = T 2 (0). Of course that it was a necessary check of consistency when 
doing the analysis without the scaling law at q 2 = 0. But, if we first extrapolate to the 
B meson mass with fl2.1p , the obtained results are just several points concentrated in the 
vicinity of the zero recoil point (q^ax = { m B — f^K*) 2 )- These values are important for 
studying these form factors in B — > K*£ + £~, but for B — > K*^ another extrapolation {i.e. 
to q 2 = 0) is needed bringing uncontrollable uncertainties to the final result. Let us also 
note that T 13 (g 2 ) do not contribute at g^ ax and, since there is no momentum injection, 
the measure of T 2 (q 2 nax ) is very clean. 

4 APE collaboration already started the lattice study of semileptonic and radiative heavy to light tran- 
sitions by implementing the new improvement program. First results are expected by the end of this 
year 
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3 Unitarity bounds 

The behavior of form factors with the values known at several points can in fact be con- 
strained by some general principles of the theory. Unitarity bounds on the form factors 
were generated for different quantities in heavy to heavy transitions [21, |26[ , as well as 
for the form factors in heavy to light case J27j, by using the old idea of [ 28| . The lattice 
data in this analysis were included first in P5| , and then in [30]. Here we give a brief 
outline of the method applying it to our problem. 
One starts from the two points correlation function: 



(3-1) 



IW(g) = i J rf 4 xe^(0|T{^(x),£ /3 (0)}|0) 

= p M - a/3 n_(g 2 ) + p+ a/3 n + ( g 2 ) 

where the tensor current is = (l/2)sa /JiV b, and projectors: 



(3.2) 



project out positive and negative parity parts which come from 1 + and 1 intermediate 
states respectively. When contracted with q v qP , this general correlator gives: 



i J d 4 xe^(0\T{j^x), jt(0)}|0> = (q,q a - q'g^U^q 2 ) 



(3.3) 



Similarilly, when we deal with jjj = (l/2)s<j^ llJ q 1/ 'y 5 b we obtain: 



i / d 4 xe^(0|T{j», ft (0)}|0) = (q,q a - q 2 9lia )U + (q 2 



(3.4) 



As usual, the functions H±(q 2 ) can be written in the form of n-subtracted dispersion 
relations: 



n ± (g 2 



n-1 



7T J t n (t-q 2 ) ^ 



(3.5) 



Unknown subtraction constants are eliminated by taking derivatives, so that we finally 
have: 



1 <9 n n ± (g 2 



ImU ± (t) 



dt. 



n! d(q 2 ) n Ti J (t-g 2 )™ +1 
The spectral functions can be obtained from the unitarity relation: 



(3.6) 



{.QfiQa - q 2 g lia )ImYl±(t 



it 



yE / ^r(2vr) 4 ( 5( P -p r )(0|j / ,(0)|r)(r| J t(0)|0), (3.7) 



where the sum runs over all possible hadronic states with suitable quantum numbers and 
with the phase space integration for each allowed state. In this sum of positive terms, the 
contribution from \BK*) is: 

/mn - (t) -T& A3/2(t) ^" t+)|Tl(t)|2 

where we have used the fact that the matrix element for BK* production is described by 
the same set of form factors but real in the region t > t+ (t± = (mg ±m^.) 2 and as usual, 
\(t) — (t — t+)(t — t-)). Note that we do not take T 3 (t) into account as mentioned at 
the beginning. Now, we may invoke duality and calculate functions xt(Q 2 ) (Q 2 = 
perturbatively. For this to be valid, we should be far away from the region where the 
current can create resonances. Assuming in this case that Q 2 = is far enough, we obtain 
to leading order: 

N c r 1 x 2 (l- x) 2 (2u-3(l- x) -3u 2 x) 
X2 { } ~ 2! \Qn 2 m 2 Jo (1 - x + xu 2 ) 2 



+ _ N c ri xHl-xn3x(u 2 -l) + 3-8u) 
Xs W " 3! 32ttX-/o (l-x + xu 2 Y dX 

where u = m s /m h and we take it to be zero or 1/25. Numerical values are then: 

( u = 0) X 2 = 2-375 x l(T 3 /ra& 

Xt = 2.375 x 10- 4 /m^ 

(u = 1) X 2 = 2-121 x 10~ 3 /m 2 b 

Xt = 1-167 x 10~ 4 /m 4 b (3.10) 

Putting this altogether in dispersion relations we arrive to the following inequalities: 

A 3 / 2 (t) 



1 r°° X^ z (t . , l2 



(m 2 R — m 2 ' 2 



B HI'K*) \ /V \ L ) I rji /,\|2 



12n 2 xt Jt + t 6 



A x / 2 (t) 



T 2 (t)| 2 < 1 (3.11) 



To extract the information on the form factors in the physically interresting region for 
B — > fT*7 decay accessible on the lattice, one performs a conformal mapping: 



1 + z = (m B + m K *) 2 -t 
1 - z V 4Nm B m K * 

by which the regions t- < t < t+ and < t < t- are mapped into the segments of the real 
■1 and z max > z > z min respectively, while two branches of the root are 
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mapped to upper and lower semicircles of \z\ 



!■ Zmin and z max 



are given by 



m B + — 2y/Nm B m 7T 
m B + + 2y/Nm B m 1T 



N-l 
yfN+l. 



where N is a free parameter whose role we explain below. It is worth noting that for 



■N 1, z m i n 



1 r 2n 



2tt Jo 



0. Generically rewritten, the transformed inequalities (|3.11|) are: 

If dz , 



dOwi^Ti^Y < 1 
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=i z 



(3.13) 



The so called outer functions Wi(8) are known and analytic on the circle, and the analytic 
functions 4>i(9) can be found as solutions of the Dirichlet's problem with |0j($)| 2 = u>i(i?) 
on the circle |31|. Their explicit forms in our case read: 
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(3.14) 



The space of analytic functions on the unit disk with 



< oo 



(3.15) 



on the boundary is the well known H 2 space. To make use of the rich theory of H p 
spaces we have to ensure our form factors to be analytic on the disk i.e. below the 



threshold for the BK* production. The problem with the poles at B* (5.42 GeV |35|]) and 
B s i (5.67 GeV |34fl) is solved by multiplying the form factors by Blaschke functions which 
simply remove the poles sending them on the circle (moduli of these functions are equal to 
one): 



T^t) -> Biz^T^t) 



Z Zpole > Ti(t) 

1 ZZpolei 



(3.16) 



Besides, the subthreshold singularities must be taken into account. This cannot be done 
in a model independent way. Some reasonable model for the cut must be employed. Such 



analyses were performed in |^5| demonstrating that the effect of these singularities is small 
and that the final bounds are just few percents weaker. While this is very important for 
this sort of analyses for some subtle quantities, like the slope and the curvature of the 
Isgur-Wise function in B — > D*iv (for the latest results in this field, see detailed analysis 
in P3|), in our case few percents are really unimportant. Now, to generate the bounds 



we may use different approaches. Here we adopt one from |H| |32|] which is particularly 
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elegant. The solution is presented in appendix A. The resulting bounds for the case when 
the form factor value is known at one point z = is: 

i (-0i(O)T 4 (O) + 0i(z)T,(*)(l - z 2 )) 2 + ^(0) 2 T,(z) 2 (l - z 2 ) < 1 (3.17) 
When the value is known at two points z = {0, Zq} the bounds are solutions of: 

, 1 x 2 |-T2 fe(0)T,(0)(2; - 2) + faizoMzoW - ^)(1 - 2 2 ) 

[Z — Zq) [ZqZ l 

-0,(^)^(^)^0(1 - z z){l - z 2 )] 2 + UzafTtizofil - z z) 2 {l - z 2 ) 
+<j )l (z) 2 T i (z) 2 (l - z z) 2 (l - z 2 ) - 2<t> i (z )<t> i {z)T i {z )T i (z){l - z 2 )(l - z 2 )(l - z z)} < 1 

(3.18) 

Now when we have the expressions for the bounds, we can apply them to the case at hand. 
To do this, let us remark that the coefficient N in ( 3.1 2|) is free and for commodity we can 



always choose it in such a way that the value of the form factor is known at z = 0. So the 
formulae for the bounds are always applicable by a simple adjustment of N. Here we use 
APE results for T\p(f) (technical details can be found in |[L3|| ). After the extrapolation to 
rriB, the results that we use are: 

T 2 (19.2) = 0.25 ±0.02 (0.26 ± 0.02) 
T 2 (16.8) = 0.24 ±0.04 (0.27 ±0.05) 

Ti(16.8) = 0.56 ±0.08 (0.59 ±0.13) (3.19) 

The above values are obtained from linear (quadratic) fit. The possibility to extrapolate 
the results obtained at nonzero recoil point (but in its vicinity) by the heavy quark scaling 
law was explained in JTJ . 



The typical bounds for the form factors T\(t) and T 2 (t) are depicted in Fig.l. The 
bounds generated by using only the lattice results in this fashion, do not give satisfactory 
results, but they can be more and more useful with the results of forthcoming simulations, 
i.e. with a bigger range of q 2 where the form factors are measured with a reasonable 
precision. For the coherent and immediate usefulness, the bounds should be formulated 
for lattice accessible masses. Then, they could even serve to interpolate to q 2 = 0. We plan 
to return to this question soon. For the moment, the bounds on the form factors in the large 
recoil region are so wide that not even a hint on the form factors' behavior can be obtained. 
However, it should be emphasized that the situation from Fig.l is expected. The g 2 -region 
accessible from this decay is so large, that one cannot expect a miracle when constraining 
by one or two points at one end of this region. Since we want to learn something more on 
the functional dependence of our form factors, we obviously need some information in the 
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Figure 1: TTie lattice constrained unitarity bounds on Ti(q 2 ) and T 2 (q 2 ) form factors. The 
APE lattice point(s) used as constraints are marked too. 
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large recoil region. Before turning to that point, let us give some more remarks. Apart 
from useful scaling laws, the symmetry of heavy flavors opened the possibility to link this 
decay to the decays generated by V — A current [PEfl - Namely, by using the equation of 
motion for the heavy quark in its rest frame (7°6 = b), we have saoib = — is / ~fi(l — 75)6. This 
allows for several relations among different form factors, valid in the small recoil region 
(the one accessible on the lattice). This fact was not used enough in the previous studies. 
It is also interesting to note that in the same limit, the authors of Ref.|37j demonstrated 
that the axial current couples only to il + and not to |1 + , which proves our assumption 
on the pole for T 2 (t) in this section. 



4 Light cone sum rules 

In seeking to answer the question on the form factor behavior, one must try to have some 
information in the large recoil region. This is the region where QCD sum rules proved 
to be successful. Here we recalculate Ti(q 2 ) and T 2 (q 2 ) paying some more attention to 
the numerics. In what follows we point out the differences with the papers where these 
calculations were originally performed. 

Ti(0) was extensively studied in the sum rule approach [^3[ |4(| [4l|]. In the case 



of B — > K*£ + £ these form factors were calculated by means of the classical sum rule 
technique [Q, as well as in the light cone (LC) approach [fj3| . The authors of Ref. ||21 



calculated the form factors for B — > plv decay and investigated the discrepancies between 
two approaches, concluding that the LC sum rules are more reliable in the large recoil 
region. LC sum rules seem to be more useful for us also because they have well defined 



heavy quark limit [23|. Brief and (in our opinion) the nicest review on LC sum rules in 
heavy to light processes, can be found in [fHJ]. In the interresting region, the recoiling s 
quark is very energetic and two competing mechanisms are on the stage: hard and soft. 
The hard one comes from the configuration where the constituents are at small transverse 
distances so that the hard gluons are exchanged. This leads to the counting rule which 
in the heavy quark limit gives Tj(0) ~ m/ 2 . The soft part comes from the so called 
end-point configuration, where almost all the hadron momentum is carried by s quark 
so that the transverse separations are not constrained any more and the nonperturbative 
technique must be used. In [j2U[| , it was shown that the soft contribution has the same 
scaling behavior (i.e. in our case Tj(0) ~ m H 3 ^ 2 ). 
The sum rule can be derived by considering the correlator: 

F^q,k) = i J rf 4 xe^(K*(fc,e A )|T{^(a;) ( T^g^(a;),6(0^75rf(0)}|0) 

= £fiiyaf3e* x u q a k l3 F(q 2 , (k + q) 2 ) (4.1) 

Formally, the hadronic representation of this correlator can be obtained by inserting the 
complete set of states. With (B(k + g)|6ry5<i|0) = fB m B/ m b an d (|1.3[), we can write down 
the dispersion relation for the amplitude F(q 2 , (k + q) 2 ): 

F(q 2 , (k + qf) = r p (£> 8 ) ds + subtraction (4.2) 
Jml s - k + a) 2 
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where the spectral density is written as the B ground state contribution plus a contribution 
coming from the higher excitations in that channel: 

p(s, q 2 ) = 2T x tf)£^5{s - m|) + p e *V, s) (4.3) 

In the same manner, by using the current (l/2)sa IJ , I/ q L 'j5b, for the second form factor we 
have: 

F^q,k) = e* {x)fl F(q 2 ,(k + q) 2 ) and 

p(s, q 2 ) = iT 2 {q 2 ){m 2 B - m ^)— - m%) + p exc (q 2 , s) (4.4) 

m b 

On the other hand, the amplitude F(F) obtained from QCD, satisfy the same dispersion 
relation (|4.2|) , starting from m 2 . Then, one invokes duality of the two representations 
in the region above a threshold Sq (which is not a priori known), and equalizes both 
representations to end up with: 

2T l (q 2 ) f B m 2 B _ 1 fo ImF QCD {q 2 , s) 



m 2 B — (k + q) 2 m b ir J m \ s — (k + q) 2 

(m B - m K ,) = - / — rfg (4.5) 



m 2 B - (k + q) 2 m b K B K tx J m 2 b s - (k + q) 

The correlator is defined in such a way that only single dispersion relation is to be used. 
In the classical approach, both mesons are treated symmetrically in a double dispersion 
relation, so that two threshold parameters are in the game. Here we have only one, and 
this can be viewed as one of the important advantages of the LC approach. QCD part can 
be calculated by expanding the T product of the currents near the light cone x 2 = for 
the large space-like momenta (k + q) 2 < (when the virtuality of the b-quark is large). 
Thus, one uses the perturbative expansion of the b-quark propagator in the external field 
of slowly varying fluctuations inside the K* meson. The remaining matrix elements define 
the light cone wave functions. To leading order (bare 6-quark propagator) one deals with 
the following wave functions [45, p3| : 

(K*(k,XM0)^d(x)\0) = -ifkM{4 X)k --< {X) K) t due lukx Mu^), (4-6) 



o 



(K*(k,X)\s(0)lA^) = K^T-fK^K* J o due lukx ct>\\{u^) 



+ - fK*m K * [ du e iuk *g? (u, (4.7) 



(K*{k, AM0) 7m 75^)|0> = ^ aP e< x »k a xf i f K .m K . £ due mkx gf '(«, /x) . (4. 



10 



u is the momentum fraction carried by the valence s quark in the K* . In the above def- 
initions, the gauge phase factor between the quark fields is implicit to ensure the gauge 
invariance. In the Fock-Schwinger gauge, above definitions are of course, exact. All dis- 
tribution amplitudes are normalized to one. Here to leading order means to leading twist 
(defined as a difference between the dimension and the Lorentz spin of an operator). To 
leading twist - 2, the distribution functions are listed in appendix B as well as numerical 
parameters taken for the analysis. It is very important to note that one of distribution 
functions that we use is qualitatively different from the corresponding one used in [^, f|3| . 
Namely, in the recent study p2| , the distribution amplitude for the transversely polarized 
p meson was shown to be broader then the original findings of argued. We take the 



functions from [22] and account for the SU(3) breaking. In Fig.2, we display the transverse 
wave function which differs from the previous studies. Using the above definitions, for the 
QCD part we obtain: 



F QCD^ q 2 



{k + qf 



f K *m K * du 



2 Jo mf — uq 2 + uum K * — u{k + q) 2 



x 



K" 



f K *m K * 



4>±(u,li) + $||0,/i) + ug { l\u,fi) + ^g [ ±{u,fi) 



.(«)/ 



and 



F QCD {q 2 ,{k + qf 



+ 



1 r 1 



9 9 9 9 

m t + 9 u m K * 



A(ml — uq 2 + uum 2 K * — u{k + q) 2 ) 

fK*m,K* du 



2 Jo ml — uq 2 + uum 2 K , — u(k + q)' 



x 



(4.9) 



(q 2 + ukq) g { l ] (u, Li) 



u \mr K ,q — {kq) A 



2{m\ — uq 2 + uum 2 K , — u(k + q) 2 ) 



+ (kq)^(u,n) + mb /«* {kq) Mu,n) 



(4.10) 



Putting this into dispersion relations and taking the Borel transform of these relations, we 
obtain the following expressions for our two form factors: 



rj, i 2\ m b f K ,m K * 

T l(2 ) = 



Am 2 B f, 



exp 



1 1 

m B — mi 



M 2 



^ du , tii i o o 9\t 

— ex Pt^7^W -m b -um K ,)\ x 
o u uM z 



ug { l\u,ij) + $n (u,li) + ^g±'(u,Li) + 



(*), 



mbf, 



f K *m K 



-(f>±(u,LM) 



e[c(u,s )} + 



/ 2 , 2 2 2 \ (a) / \ 

-{q +m b -u m K *)gy{u,Li) 



uM 2 



Q[c(u, s )] + S[c(u, s )} 



(4.11) 



11 



0.0 



0.2 



0.5 



0.8 



1.0 



0.0 



0.2 



0.5 



0.8 



1.0 



Figure 2: Distribution amplitude <p±{u) (bold curves) is plotted for the case fi = lGeV 
(left), and fi = fif, (right). In the right figure, the two other solid curves show the result of 
the variation of the parameters a^siHb) listed in App.B. Dashed curve: in the left figure, 
it shows this distribution amplitude for the case of p meson as obtained in [2~2~|; in the right 



figure, it is 4>±(u) used in 
6w(l — u). 



Dotted curve in both figures shows the asymptotic shape 



m b f K *m K * 



exp 



m\ 



M 2 



Jo ^ 6XP W (g -™b-vm K *)} 



x 



[m h 



m K ,u 



(ug { l\u,fi) 



<$>\\(u,fi)) + 



m b f 



f K *m K * 



fi)) + 2uq 2 g { l\u,fi) 



0[c(u,s o )] + 



( 2 2 

{m b - q 



u m 



K" 



2 2 2 

u m K ,q 



9± (u,fi) 



uM 2 



Q[c(u, s )] + 5[c(m, s )] 



where M 2 is the Borel parameter, u = 1 — u and c(w, so) = wso — m 2 + g 2 w — uum 2 K , . Note 
that we used the continuum subtraction prescription explained in |^TJ. Apart from that, 

m 



the above formulae agree with 



The resulting form factors are shown in Fig 3. 



The results of the sum rule analysis are expected to be reliable for q — mi negative 



enough. For the case at hand, the range of reliability is below q 2 ~ 15 GeV 2 f23| 
From (4.11) and (4.12) (see also Fig.3), we extract the value of the form factor: 



T(0) = 0.15 ± 0.01 ± 0.02 ± 0.02 



tii 



(4.13) 



where the first error comes from the variation of m b and the sum rule parameters (Borel 
and the threshold parameters), the second error is related to the variation of the wave 
functions' coefficients (see App.B) and the third is 'theoretical' uncertainty of order 15% 
(not included higher twist and radiative corrections). To be on a safe side, we varied all 
parameters that enter the calculation in their largest reasonable ranges. This is particularly 
the case for a^~(/i) coefficient whose variation mainly contributes to the second error. Note 



that on Fig. 2 we displayed 



for only one ('central') choice of (1 GeV). As it can be 



noticed, although with different shapes of wave function, numerical consequences on the 



(4.12) 
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Figure 3: LC sum rule prediction for the Ti(q 2 ) and T 2 (q 2 ) form factor. Dashed lines 
represent the variation of the sum rule parameters as indicated in appendix B. 
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form factor T(0), when compared to are of rather minor importance (see Tab.l). This 
is due to the fact that the value of f^, that we use is smaller than the one used in [^, £43 



which compensates the higher value of (see App.B). Let us also give some cross results. 



Without SU(3) breaking, i.e. with the wave functions from [22] the above sum rules yield 



th 



T(0) = 0.13 ± 0.01 ± 0.02 ± 0.02 
We also evaluated the sum rule for B — > ptQ and the result is: 

T(0) = 0.12 ± 0.01 ± 0.01 ± 0.02 th 



(4.14) 



(4.15) 



Of course the true result for this decay rate can be obtained only when long distance effects 
from effective hamiltonian are included p8 |. In the limit where they are neglected, we have 



B{B -> fry) 
B(B -> K*~i) 



Vt 



td 



t.s 



r 2 (0) 



m b + m d f m 2 B -m 2 p \ 3 



mt + mi 



m 



B 



m 



(4.16) 



K" 



where for r(0) = \T b ^p(0)\/\T b ^ k * {0)\ we obtain r(0) = 0.776 ± 0.073. 



5 Bounds with additional constraint 

In the previous section we obtained the form factor at q 2 = 0, the value of which is 
estimated by including various uncertainties. Now we can generate the bounds from sec. 
3 with one lattice result at q 2 ~ q^ax an d LC sum rule result at q 2 = 0, to see if something 
more can be learned about the form factor behavior. These bounds are depicted in Fig 4. 
Now, the bounds are much better. It is obvious that for a clear conclusion on the form 
factors' behavior, we need some more points to constrain more the intermediate region. For 
this, the lattice results are needed if we want to treat both methods on equal footing. There 
are other sources of improvements of this analysis: for the functions x^) corresponding two 
loop expressions can be calculated; try to seriously examine the effects of subthreshold 
singularities. At the same time, the bounds need to be relaxed for the value of systematic 
uncertainties. Some of these uncertainties for the lattice results are included using the 
quadratic instead of the linear fit to extrapolate to JTjj]. 

However, even now from the two available results, it is obvious that the fully flat 
behavior of T 2 can be discarded. This does not contradict [ I3| , where the best fit was the 



pole one, with the pole mass as a free parameter. It is unlike that the nearest pole behavior 
can work in the full (large) region for this decay. This would lead to the oversimplification 
of the contributions of all possible singularities that lay behind the nearest pole, especially 
in the large recoil region. We do not know strengths nor the signs of these singularities, but 
the least we can do is to have doubts on their mutual cancellations and leave "m po i e " as a 
free parameter which is to be fixed in fit with the data. In the light of that remark, we note 

5 In H], first candidates for B — > pj have been extracted. 
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Figure 4: Unitarity bounds for the form factors Ti(q 2 ) and T 2 (q 2 ) constrained by the lattice 
result at high q 2 and by the LC sum rule result at q 2 = 0. 
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Figure 5: Form factor behavior: Bold solid lines show the double pole (pole) behavior of 
Ti(q 2 ) (T 2 (q 2 )) with T(0) = 0.15, m< Aivole > = 5.8 GeV, m< po i e > = 6.7 GeV. Dot-dashed bold 
line shows the nearest pole dominance for Ti(q 2 ) compatible with ( |2.1| ) andT 2 (q 2 ) being flat 
(dotted lines). Unitarity bounds forT\(q 2 ) are represented by dashed lines. 



that the pole fit with m po i e ~ 7 GeV reported in [T4[ is compatible with our bounds. In 
fact the dipole/pole behavior is favored (of course not with the nearest resonance masses). 
Specifically, in Figla. we can make the pole fit with Ti(0) > 0.23 and m po i e ~ mg. 
while with the result from Fig4a. this cannot be done and the double pole fit works. 
For instance, with 71^(0) = 0.15, the free parameter m^ipoie' = (5.6 — 6.2) GeV and 
m-'poie' = (6.4 — 7.0) GeV, where the variations of free parameters allow to stay within 
lattice constrained bounds in the other end of the q 2 region. This is illustrated in Fig. 5. 
That is what we obtain concerning the q 2 behavior of our form factors which is important 
for the study of the B — > K*£ + £~ (i.e. for the off-shell photon). On the other hand, the 
extrapolation of the lattice results in mass at q 2 = is surely a better way to extract the 
value of the form factor for B — > K*j. As we already mentioned such a scaling law is 
now available. Here, we will not go through the lattice analysis again, but rather make an 
important consistency check of the two approaches that we use in this paper. First, one 
can calculate the heavy quark limit of the sum rules p3fl . We recalculated that too, and 
for very generous variations of the sum rule parameters we obtain: 



T(0)m^ /2 = (2.7 ±0.3) GeV 3/2 . 



(5.1) 
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Now we can make the fit with the available lattice dataf] as: 



i% 2 = (2.7 ±0.3) GeV 3/2 (l 
for lattice attainable m#, extrapolated to q 2 = (with m~ 3 / 2 scaling law, see table2. in 



T(0)m^ z = (2.7 ±0.3) GeV' A ' 2 [1 + — + JL] ( 5 . 2 ) 



13|]) to fix A and B. The values of the coefficients obtained in this way are: 

A = -(2.1 ± 0.2) GeV B = (1.0 ± 0.2) GeV 2 . (5.3) 

These numbers are big, which indicates that the corrections to the infinite quark mass limit 
are large. Now, we simply extrapolate to the physical B mass and the resulting value of 
the form factor is: 



T(Q,m B ) = 0.14 ±0.02. (5.4) 

If we include ~ 10% of errors to this result due to quenching, we obtain the lower bound 
as T(0, m#) > 0.11. So we can reduce the uncertainty on T(0) reported in the previous 
section ending up with the still conservative result: 

T(0) = 0.15 ±0.04. (5.5) 

Another consistency check of LC sum rules and lattice predictions can be made for the 
strength of the nearest pole contribution, i.e. the residue ~ fg, Qbibk* which supposedly 
dominate the form factor Ti(q 2 ) behavior in the zero recoil region. We did not calculate 
that but plan to do it with the new lattice data. 

In the tablel., we enumerated some of the existing results insisting on the QCD directly 
based predictions. For completeness (comparison), some quark model predictions are also 
given. 

6 where we again use APE lattice results 
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Reference 



T(0) 



2 > 
max i 



This work 



0.15 ±0.04 0.26 ±0.02 13 



BHS (H 
UKQCD H 
APE PI 



LANL 15 



0.10 ±0.03 

o.i5l:K 

0.09 ±0.02 
0.09 ±0.01 



0.33 ±0.07 

0.271$ 
0.26 ±0.02 
0.23 ±0.01 



Colangelo a/.[|2| 
Aliev et al. |4j 



Ball |E 

abs ; 

Narison 



23 



10 



DPR 41 






19 ±0.03 


0.14 ±0.02 




0.15 







19 ±0.03 







16 ±0.03 







15 ±0.02 






0.20 





MNS U 
Stech m 



Olsson, Veseli 



0.18 ±0.05 

0.18 
0.13 ±0.03 



0.25-0.35 
0.35 



Table 1: The form factor T(0) and T 2 {q^ nax ) obtained on the lattice (quoted lattice results 
are obtained with pole/dipole scaling law), from LC and traditional sum rules, and quark 
models. 



6 Conclusion 

In conclusion, let us briefly outline the result of our analysis. We generated the unitarity 
bounds for the B —>■ K*^ form factors accessible on the lattice. The bounds constrained 
by lattice results only are not very restrictive. They will become more useful with the new 
lattice data. Besides the heavy quark scaling laws in the small recoil region, the scaling 
law at maximum recoil is accounted for. The former contradicts both form factors being 
pole dominated because of the condition Ti(0) = ^(O). To have some more insight in the 
form factors' behavior that can be drawn from this analysis, we also recalculated the light 
cone sum rules for them. By including the errors of theoretical uncertainties and varying 
the sum rule parameters in their large reasonable ranges, we obtain that pole/double pole 
fit is in fact better. For the consistency check, we also evaluated the sum rule in the 
heavy quark limit and fixed parameters for the extrapolation to the B-meson mass by the 
existing lattice results. Such obtained results are very consistent with the results (bounds) 
obtained from LC sum rule. In generating the bounds, the pole that exists for each form 
factor is properly accounted for. Then, the result that one of the form factors may be fitted 
by dipole may seem a bit surprising. This is not worrisome at all since the unconstrained 
bounds are very wide (for instance |T\ (0) | < 10), and the decisive limitations come actually 
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from the results taken to constrain these bounds. At the end, the result that we extract 
from this analysis for the form factors at q 2 = and with all constraints included, falls in 
the conservative range 0.15 ± 0.04. With this, we can predict the hadronization ratio: 



R 



K* 



(12.2 ±6.6)% 



tub m 2 B ~ m \* 



m b 



m 2 B 



IT 



(6.1) 



with the above formula being valid to leading order in QCD |48|]. We also give our prediction 
for the branching ratio for this decay, with C 7 = —0.333 and r B = 1.6 ps [p5| : 



B(B -> K*j) = \V*V tb \ 2 (2.41 ± 1.27) x 10" 2 
= (4.5 ±3.2) x 10" 5 



(6.2) 



|V^*V t ft| is also taken from ]3"5| . With these parameters, we can extract the experimental 
values of the form factor: T(0) = 0.17±0.06. New theoretical and computing improvements 
are thus more than needed. 
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Appendices 



A Explicit derivation of the unitarity bounds 

In this appendix we outline the way to get the unitarity bounds on the form factors with 
their values known at several points as constraints. The 'pure' unitarity bounds can be 



obtained as usual (see any of references ^5], [28], |29|, |30|): 

\T(z)\ < ' ' 



(*)|VT 



(A.l) 



and they are very weak for the case we study. Suppose now that form factor is known 
at several points {z{}i . Recall that the functions we consider are real and analytic in 



1 belonging to unit sphere of TCI 



i.e.: 



II Tw \\ L 2< 1 (A.2) 
For example, when the form factor is known at z — {0, Zq}, we look for the analytic function 



which takes the two prescribed values and satisfies the constraint |A.2j , which is equivalent 
to methods used in ffijfl . For this to work, the sequence of points Zi should be uniformly 
separated [[31]]. If this is not fully satisfied, the probabilistic approach can be adopted as it 
was done in Ref . PS| . 

The most general functional form which satisfies above mentioned constraints is: 



T{z) = ±T{ Zl )^- + B(z)g(z) 
which can be simply checked employing the usual definitions: 

n 

B(z) = J[f B,[: 



" z — Zi „ , s -X z — Zi 



i=l 



ZZi 



ZZ; 



(A.3) 



(A.4) 



Obviously, the first term in T(z) satisfies both constraints (i.e. prescribed values and ( |A.2|) 
) while the second term contains an unknown function g(z). The aim is then to look for 
the g(z) belonging to H^, closest to the given kernel i.e. |32|, [26| : 



where the kernel: 



V 1 = (min ||^-^||) ; 



. " T(zi) 1 - 



r{ Bi(zi) z- zi 



(A.5) 



(A.6) 



is just the first term of T(z) divided by B(z). It is convenient to solve this problem by 
considering its dual [|3~T|: 



/i 2 = sup |— r dQ G(6)il>(6)<f>{6)\ 

Ge-H" 2 ,||G||<1 27T Jo 



(A.7) 
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G{z) = GkZ k is surely inside the unit sphere of H 2 if G\ < 1 which again reflects 
the constraint (|A.2|) . Then we simply apply the residua theorem to calculate integrals and 
obtain: 

oo 

// = sup £^d n) (A.8) 

where by £(n) k we compactly wrote the result of integration of each term of the power 
series, and n stands for the general case of n — 1 known form factor values at fixed z = z n 
which we then vary to obtain the functional dependence that we are interested in. Now, 
using the Schwartz inequality, we eventually obtain: 



fc=0 



1/2 

< 1 (A.9) 



It is very easy to show that for n = 1 we obtain the pure unitarity bounds ( |A.1| ). 

For example, when n = 3 i.e. the form factor value is known at Z\ % % = {0, zq} and we 
look for its functional dependence at z% = z, in ( A.8 ) we have: 




+ Tjzp) 1 - 2021 -^ + T{z) l - Z « Zl - z2 \ G 

ZqZ Zq — Z Zq Z — Zq Z I 

z z)(l - zl)z k - 1 + -^-(1 - z z)(l - z 2 )z k A G k (A.10) 

Z — Zq I 



and by inserting it to (|A.9|) , we obtain the bounds given in Eq.(EyJ). The case n = 2 that 
appears in the paper, is obtained in very much the same way. 

Note that each time the form factor is multiplied by a proper Blaschke factor to account 
for the pole. In order not to make it heavier then necessary, this is implicit in all the 
formulae displayed so far. To have some feeling about the numerical analysis, we give some 
significant numbers: 

• For the form factor Ti(q 2 ), the measured lattice point is at q 2 = 16.8 GeV 2 — > z = 
when = 1.130 and the pole position in z-plane is at q 2 = 29.38 GeV 2 — > z po i e i = 
-0.2188. Also z(q 2 = 0) = 0.1443 was used. 

• For the form factor T2(q 2 ), the measured lattice point is at q^ax = 19-2 GeV 2 — > 
z = with N — 1; the pole position q 2 = 32.15 GeV 2 -> z pole2 = -0.2811 and 
z(q 2 = 0) = 0.1742. 
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B Distribution amplitudes and numerical parameters 
used in LC sum rules 



For completeness, we display the wave functions that we use in the evaluation of the sum 
rules. We show only the steps which are particularly important for our analysis. The rest 
can be found in e.g. ||22j| . To leading twist accuracy, there are two functions that are 
independent: 



^l,|| ( M ) A*) = 6tt(l — u) 



.(B.l) 



To the same accuracy, the other functions that figure in Eq.(4.6,4.7,4.8), i.e. g±' a \&\\ are 
not independent. They are expressed in terms of <p\\(u, p) (see f22] , |2^|). Their explicit 
forms read as follows: 



9 ( ?\u, p) = \ (l + C 2 ) + + 4(a0^ (15C 4 - 6C 2 - l) + 4(p) 



en 



2 V5 



9±\u,p) = 6u(l - u) 



C + 4(/x)l (2C 2 - l) + a«(A*)| (C 2 - + gS(aO^ (f C 4 - 9( 2 + 1 

(B.2 



where £ = 2-u — 1. In the case of the p meson ai = = 0, and a^"(l GeV) = 1.50 ± 0.75, 



a 2 (l GeV) = 1.35 ± 0.75 p2 |. The coefficients ai^ give rise to antisymmetric parts of the 
wave functions which acquire nonzero values upon SU(3) breaking. To obtain the values 



of the coefficients, we use the results of |^5| for corresponding moments (C n ) : 

K" 



(C 2 ); 



0.75, 



(C 2 )i 



0.85, 



(B.3) 



for a 2 {lGeV), and 



(C 3 



(0 



j-.ii 



0.50 - 0.55, (C)^ 11 = 0.15-0.20, 



;b.4) 



for 0^3(1 GeV). Thus obtained values are then evoluted to pi ~ 5 GeV, so that we obtain 
the following results (ranges of values) for the coefficients: 

a\(p b ) = 0.64 ± 0.09, al(p b ) = 0.31 ± 0.37, a|(/x 6 ) = 0.38 ± 0.15; 
a^(p b ) = 0.69 ±0.10, 40&) = 0.21 ± 0.41, a^(p b ) = 0.42 ± 0.16. (B.5) 
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The values of these coefficients are varied in their respective ranges and the form factor 
results are relatively insensitive to it, as we reported in (4.13). This is true provided 
SU(3) breaking is carried also for the constants fy and especially fy. Again, from p2f : 
f e p xp = (205 ± 10) MeV, ft = (160 ± 10) MeV, and with Uk* / f^? = 11 fr o m [@> we 
have: 

f K , = (215 ± 10) MeV ft, = (170 ± 10) MeV (B.6) 

So far the coefficients cij and ft (V stands for a vector meson) were obtained in the QCD 
sum rule approach. There are very few lattice calculations which were performed using 
small lattices and a relatively poor statistics. For the case of pion distribution amplitude, 
the authors of Ref.[^9| studied the second moment for which we can extract the value 



of a^ilGeV) = 0.44 ± 0.53. The errors are very large and they are mostly due to the 
modest statistics. Only recently this result was improved in the study of ||50|in from which 
we extract 0^(1 GeV) = 0.40 ± 0.27. It is rather difficult to compare these results to that 
obtained from the QCD sum rules a£(l GeV) = 0.43, although they agree well at first 
sight. The sum rule study |45| of the leading twist distribution amplitude for the pion 



needs also to be updated as it was done for the case of the rho-meson |22[ . 
In numerical analysis of the sum rules, the following values are used: 

• b quark mass varied - m h = (4.7 ± 0.1) GeV 

• threshold parameter varied - Sq = (34 ± 1) GeV 

• Borel parameter M 2 = (6 — 10) GeV 2 

For the sum rule to estimate the value of fs-, also (qq)(l GeV) = (—245 ± lOMeV) 3 , 
and {qgaGq){l GeV) = 0.65 GeV 2 (qq)(l GeV) was used (the values taken from |2l||). 
The sum rule for fs, is taken from |SJ where the radiative corrections are not included. 
This hopefully cancels some uncertainties since these corrections (though needed) are not 
calculated for the form factors yetQ. The stability of the sum rules when the Borel parameter 



is varied is shown in Fig. 6. Also the change due to the subtraction plU , with and without 
delta-function term is shown on this figure. 



7 In this reference the authors also calculated the moments for the p meson. The lattice estimate of fy 
and the comparison of these results with the results of QCD sum rules will be addressed elsewhere. 

8 While this paper was in writting, the authors of Ref.|5l| calculated the radiative corrections to the LC 
sum rule for the f + (q 2 ) form factor which enters B — > ttIv. Indeed, their results show that the uncertainties 
due to these corrections almost exactly cancel against the corresponding corrections for fs- 
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Figure 6: In the figures a) and b) the form factors' dependence on the Borel parameter is 
plotted. Dashed lines show the variations of and sq. In the figures c) and d) we plot 
the form factor behavior for Ti(q 2 ) and T 2 (q 2 ) for M 2 = 8 GeV 2 and rrib = 4.8 GeV with 
and without (dashed) delta function term for the continuum subtraction in LC sum rules. 
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